#### Figure S10: Bivariate Support for Censorship ####

# Libraries
# library(here)
# library(rio)
# library(tidyverse)
# library(ggthemes)
# library(lubridate)
# library(estimatr)

# data_pnas = import(here("Data","data_pnas.rds"))
# source(here("Code","Functions","funcs_theme.R"))
# source(here("Code","Functions","funcs_bvplot.R"))

fig_s10 = plot_bv_maga(data_pnas, norm_censorshipc)

print(fig_s10)

ggsave(here("Plots","Supplementary","figure_s10.png"),
       fig_s10,
       units = "in",
       width = 6, height = 6,
       dpi = 600)